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It  is  well  known  (see  Stratonovitch,  1963,  for  instance)  that  one  can  as¬ 
sociate,  under  some  conditions,  to  a  stochastic  continuous  automata  (i.e..  a 
stochastic  differential  equation)  a  so-called  Fokker-Planck  (F-P)  equation  in 
the  probability  distribution  of  the  state  variables.  In  this  note,  we  wish  to 
characterize  conditions  under  which  the  F-P  equation  admits  a  stationary 
solution  of  the  Gibbs  type. 

Let  x  a  n-dimensional  vector  of  state  variables,  and  W(x,t)  the  proba¬ 
bility  distribution  of  the  state  variables  described  by  x  at  time  t.  The  F-P 
equation  is: 
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where  da(x)  is  the  drift  vector  and  Kals(x)  is  the  diffusion  matrix  ( see  Stratonovitch, 
1963,  p.  76). 

The  stationary  solution  w(x)  of  the  F-P  satisfies  the  equation: 

£  -  0  ,  (1) 

a=l  UXa 

where  we  have  defined  the  probability  current  Ga(x): 

Ga{x)  =  da(x)w(x)  -  ~if„o(x)ti)(x)  .  (2) 

l  OXq 

In  order  to  find  the  stationary  solution,  we  do  not  assume,  as  Stratonovitch 
and  everybody  else  does,  that  Ga(x)  =  0  and  set  w(x)  =  e~u<-x^  in  equation 
(1),  obtaining: 


- 1  jrfuw  (<■<*> +  -  \-kK'“w)  - 0 

Assuming  that  the  diffusion  matrix  is  constant,  that  is  Kaa{x)  =  Kal3,  we 
obtain: 
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=  -VU{ d  +  VC/)  +  V(d  +  VU)  =  0  . 

We  finally  obtain 

(V-  VC/).(W  +  d)  =  0,  (3) 

which  is  the  condition  for  stationary  distribution. 

Thus,  one  solution  is: 

VC/  -j-  d  =  0  <=  d  =  — VC7,  (4) 

which  is  equivalent  to  the  so  called  potential  conditions ,  that  amount  to  say 
that  d  is  the  gradient  of  a  potential.  If  the  potential  conditions  are  satisfied 
the  probability  current  Ga(x)  is  identically  zero,  and  thus  detailed  balance 
holds.  Therefore  we  recover  the  well  known  result  that  detailed  balance  im¬ 
plies  the  existence  of  a  stationary  Gibbs  distribution  u>(x)  =  e_cr(x).  However 
condition  (3)  shows  that  the  converse  is  not  true.  In  fact  equation  (3)  has 
also  the  solution 

( V  -  VC7)  ■  f  =  0,  (5) 

with  f  =  VC/’  +  d  and  this  solution  is  not  trivial  only  iff  ^  0,  that  is  if  d  is 
not  the  gradient  of  a  function.  Equation  (3)  has  therefore  a  ularger”  space  of 
solutions  than  the  one  represented  by  the  potential  conditions.  Of  course  in 
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both  cases  the  solution  U  must  be  such  that  tu(x)  =  e  U{x)  is  a  probability 
distribution,  and  therefore  the  following  additional  condition  must  hold: 

J  dx  e~u^  <  oo 

A  simple  and  interesting  example  that  proves  the  existence  of  non-trivial 
solutions  U  such  that  u>(x)  =  e~u ^  is  the  following. 

Example  of  existence 

Consider  the  stochastic  differential  equation  in  TV 

f  x  =  -2 x  +  y  +  (x(t) 

\y=-x~2  y  +  tv{t)  lb) 

where  £x(£)  and  £y(£)  are  Gaussian  noise  terms,  that  is 

<  >  =  <  &(*)&(*')  >  =  -  0  • 

The  F-P  equation  associated  to  (6)  is 

jjWfot)  =  -V-(U'(x,l)d(x))  +  2V1W(x,t) 

where  the  drift  vector  is  d(x)  =  (—2x  +  y,—z  —  2 y).  It  is  easy  to  verify 
that  d(x)  is  not  a  conservative  field,  so  that  detailed  balance  does  not  hold. 
However  a  stationary  solution  of  the  F-P  exists,  with  u’(x)  =  e~u^  and 

U{x)  =  xJ  -) -y2  . 

In  fact,  defining  f  =  VU  +  d  we  have 

f  =  (2x  -  (2x  -  y),2y  -  (x  +  2y))  =  (y,  -x) 
and  therefore  equation  (5)  is  satisfied,  since 

(V  -  VU)  ■  f  =  V  •  f  -  VU  ■  f  =  2 (x,y)  ■  (. y ,  -x)  =  0  . 

Notice  that  in  absence  of  noise  the  differential  equation  (6)  is  linear,  with 
characteristic  eigenvalues  A  =  —  2  ±  i,  and  the  associated  trajectories  are 
inward  spirals.  This  makes  perfectly  plausible  the  fact  that,  in  presence  of 
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noise,  the  probability  distribution  of  the  variables  is  a  Gaussian  centered  in 
the  origin. 

Remarks: 

•  In  the  linear  case,  that  is  when  d(x)  =  Ax  and  A  is  a  symmetric  matrix, 
detailed  balance  always  holds,  because  d  =  -V(7  with  U{x)  =  -  jx.4x. 
However  the  stationary  solution  u>(x)  =  e~uW  exists  only  if  the  matrix 
A  is  negative  definite,  that  is  if  u>(x)  is  integrable. 

•  Stratonovitch  (1963,  p.  79)  says  that  even  when  potential  conditions 
are  not  met  but  d  is  a  linear  function  of  x  and  Kag(x)  are  independent 
of  i,  the  F-P  equations  can  be  solved.  In  fact  it  is  easy  to  see  that 
if  d(x)  =  Ax  and  A  U  not  symmetric  the  potential  conditions  do  not 
hold  but  the  function  U(x)  =  —  IxAx  is  a  solution  of  equation  (5). 

•  If  the  forces.  da  in  the  Langevin  equation  are  conservative,  i.e.,  d  = 
—  VU,  then,  if  the  fluctuations  are  thermic-like,  detailed  balance  is  sat¬ 
isfied  and  a  Gibbs  stationary  distribution  exists  (Equation  4  is  satis¬ 
fied). 

•  It  appears  that  our  results  may  be  derivable  from  the  formulation  of 
Graham  (1980)  and  the  more  general  case  considered  by  Jauslin  (1984) 
and  Zeeman  (1988).  An  in-depth  analysis  of  many  properties  of  the 
Fokker- Planck  equation  relevant  for  this  note  can  be  found  in  Tan  and 
Wyatt  (1985). 

Acknowledgements  We  are  grateful  to  J.  Wyatt  who  read  this  note 
and  provided  several  remarks  that  go  well  beyond  the  scope  of  it. 
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